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Few-cycle soliton propagation 
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Soliton propagation is usually described in the "slowly varying envelope approximation" (SVEA) regime, 
which is not applicable for ultrashort pulses. We present theoretical results and numerical simulations for both 
NLS and parametric (x^^^) ultrashort solitons in the "generalised few-cycle envelope approximation" (GFEA) 
regime, demonstrating their altered propagation. 
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I. INTRODUCTION 

This paper demonstrates the behavior of both Kerr and 
parametric soUtons in the few-cycle regime. We compare the 
usual "many cycle" or "narrowband" Slowly Varying Enve- 
lope Approximation" (SVEA) descriptions to those obtained 
with the more accurate "Generalised Few-cycle Envelope Ap- 
proximation" (GFEA) theory. We present a range of results 
using both analytic and numerical methods. 

There have been many suggestions for using solitons in 
communications as well as optical logic (see e.g. [|l], ^ |[]). 
These generally emphasise the stability of soliton profiles over 
long propagation distances, and the potential for higher data 
rates with shorter solitons [0, |[|. It is clear that the applica- 
bility of the nonlinear Schroedinger (NLS) equation to Kerr 
soliton propagation is only valid in the SVEA regime, where 
the soliton pulse contains many optical cycles. However, since 
the shorter the pulse, the fewer optical cycles it contains, the 
drive to ultrashort solitons will eventually reach the few-cycle 
regime. 

We discuss the change in propagation properties of few- 
cycle solitons when the GFEA theory is used. A description 
of the GFEA propagation equation is presented in and a 
detailed derivation is given in [0]. It emerges that few-cycle 
pulses receive an extra "phase twist" compared with many- 
cycle pulses, and this raises the question of how the funda- 
mental characteristics of soliton propagation might be pre- 
served in the few-cycle regime. 

In simple cases, an SVEA propagation equation can be con- 
verted to the corresponding GFEA form, accurate in the few- 
cycle limit, merely by applying the following operator to the 
polarisation term - 



{l+ldr/Ojof 



coo 

,a(l-a)2 4 + 0(a>4) 



(1) 



Here 3, is used as a compact notation for d/dt; (Oo is the 
carrier frequency, and a is the ratio of the group velocity to 
the phase velocity. 



II. NLS SOLITONS 

NLS solitons are hyperbolic secant (sech shaped) pulses 
that rely on the interplay of third order x'^' (Kerr) nonlinear- 
ity and the material dispersion to propagate without chang- 
ing. It is relatively simple to derive the necessary "nonlin- 
ear Schroedinger" (NLS) equation for optical pulses in a Kerr 
nonlinear medium in the SVEA limit. As usual we write the 
field in the form E{t,z) ^ A{t,z)e-"^o' + A*(f ,z)e'"o', where 
the envelope A varies slowly in comparison to the carrier pe- 
riod. 

The lowest order (bright) soliton solution of the SVEA NLS 
equation is a hyperbolic secant pulse, which at z = will be 



A(r,0) = r| sech(r|f) 



(2) 



If instead we apply the more general GFEA to the propaga- 
tion of optical pulses in a Kerr medium, we find that the pulse 
envelope of the soliton evolves according to: 



d,A + ^id,A+'-Y^fA 



= a 



(1 +;a3,/coo) 



\A\^A (3) 



l + (2-o)i 



COo 



\A\^A.(4) 



We get the approximate form in eqn.(^ by truncating the 
expansion of the GFEA operator given in eqn.(n]) to first order 
in 3, terms. In the SVEA-equivalent limit where]8rA | /coo ^ 
limit, the RHS nonlinear term is just 1% \A\ A, and the equa- 
tion reduces to the standard NLS equation. The group ve- 
locity to phase velocity ratio is a ~ coopi/Po, where p„ = 

d"m/dwr\^. 

Readers will notice that the mathematical form of the first 
order GFEA correction to eqn.(|^) appears similar to the so- 
called "optical shock" terms already discussed in the litera- 
ture (e.g. see Zaspel[^ ^], or Potasek [^), which represents 
an intensity dependent group velocity. A similar term also 
appears in Biswas & Aceves [ p^ who describe it as the "self- 
steepening term for short pulses" (see also [ pT| , p^), which 
originates from a high-order dispersion effect, and is only in- 
directly a "short pulse" effect. The origins for these terms are 
not the same as for the GFEA correction terms, although they 
have the same self-steepening effect. 

Note the difference between our GFEA few-cycle terms and 
those given by the SEWA of Brabec and Krausz[pj|, which 
is due to their approximation of o = 1. Close inspection of 
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an alternative propagation equation valid for few-cycle pulses 
given by Trippenbach et. al. [ [l4| ] also reveals a (a — 2) cor- 
rection, present in the third nonlinear term on the RHS of their 
eqn.(14). 



A. GFEA Simulations 

The code developed for modelling the propagation of op- 
tical pulses in either the SVEA or GFEA regimes is an im- 
proved version of the one used in [9]. Normalised units were 
based on a carrier frequency of coo = 27iv = 207t/A, unit non- 
linear coefficient (%= 1), pulse width (r| = 1), and peak am- 
plitude Ao — I. In this scheme, A = 1 and A = 3 solitons 
contain about 10 cycles and 3 cycles respectively. Distances 
are normalised by the dispersion distance, so ^ = z(r|/p2) 
Our simulations implement the full GFEA correction to the 
polarization, as shown on the LHS of eqn.(|l|). 

A typical result for single soUton propagation is shown in 
fig. [l|. Whereas under the SVEA the pulse would be station- 
ary within the display, there is now a drift to larger f values 
(reduce group velocity) arising from few-cycle effects. The 
results of simulations over a range of values of A and a are 
shown in fig. 2. The presence of the (2 — a) term in eqn.(|l]) 
means that, to first order, the few-cycle correction vanishes 
near a — 2 and the sign of the drift reverses at this point. 
Moreover, as the bandwidth represented by A gets larger (and 
the number of cycles correspondingly fewer), the velocity 
change becomes more pronounced. Note also that, for o = 0, 
only two correction terms appear in eqn.(|l]). 




FIG. 1 : GFEA soliton propagation. This shows the pump field pro- 
file (i.e. |Ap) for a A = 3 pulse with group/phase velocity ratio of 
o = 1. The pulse starts with an offset of Tq = — 5. 



B. Theoretical Group Velocity Shift 

It is not necessary to rely on computer simulations to pre- 
dict this few-cycle shift in the group velocity, at least in the 
case of weak few-cycle effects. Biswas & Aceves [ [l0[ | have 
already provided a multi-scale method giving the effects of 
various perturbations on standard NLS soliton propagation. 
Their eqn.(13) applied to our NLS equation gives the velocity 
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FIG. 2: GFEA soliton velocity shifts (— 8v). as a function of band- 
width multiplier A. The pulses at A = 1 are approximately 10 cycles 
long. The symbols represent the velocity shifts obtained by numer- 
ical simulation, the solid lines the corresponding theoretical predic- 
tion. The numerical labels give the group/phase velocity ratio a cor- 
responding to the nearest data points. 



shift as 



VA 



dz /+rA*A dt ' 



(5) 



which we can evaluate by inserting the SVEA soliton profile 
from eqn.(^ into the few-cycle perturbation term in eqn.(^, 
namely 



R = X (2 -o)3,A*A Vcoo. 



(6) 



This is of the same form as the A, "self-steepening" term 
A,3( l^p^ in Biswas & Aceves |]lO|], allowing for the changed 
notation and different prefactors. Solving to first order in 
the perturbation and using the intermediate quantity T = 

S-^f^z (A*A) dt gives 



coo 2ri 



VA 



2-a 7 
11"- 

coo 



(7) 
(8) 



With our parameters, this gives a velocity shift of 



5v = va = — 



(2 -a) A 
2071: 



-0.01592(2-0) A. (9) 



Since we assume a small perturbation (eqn.(|g)), it is clear 
that this prediction is most valid for higher carrier frequen- 
cies and/or weaker nonlinearities, i.e. where the effect of the 
nonlinearity is small over the time of an optical cycle. 

The predictions of eqn.(^ are plotted as solid lines in fig. 
^ and the agreement with the numerical simulations is seen 
to be remarkably good, even for very wideband pulses (e.g. 
the ~1 cycle cases where A = 8). However, the effect of the 
higher order GFEA contributions, present in the simulations 
but not in eqn.(^, becomes visible near a = 2, where the first 
order corrections become small. 



ni. PARAMETRIC SOLITONS 

Parametric solitons, otherwise known as "simultons" or 
"quadratic solitons", rely on the interplay between dispersion 
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and a second order x'^' interaction to maintain fixed envelope 
profiles between a pair of pulses propagating in tandem. 

If we modify the standard (dimensionless) propagation 
equation of Werner & Drummond[|l5l to include few-cycle 
terms to first order, with carrier frequencies cOy = 2c0(|, and 
wavevectors poy = 2po0, we get 



3^(1) + - sgn(P2(t,)3?(l) 



l + i 3i 



CO,j, 



(10) 



(11) 



Here piy,Pii|) and P2v,P2<^ are the group velocities and 
group velocity dispersions respectively. We work in a co- 
moving frame where piy = pi,],, so the two pulses remain 
co-propagating. This ensures that their group-phase velocity 
ratios are identical (o Oy = 0(|,). Distance and times are nor- 

mahsed using ' = \x<Pq\ = |P2(^| /Iq, and ^ = z/zo, i = f/fo- 
In the SVEA limit (|3xX|/| /cOy 0, and |3x(|)| /co,|, 0), the 
standard anzatz gives the solutions 



^{z,t) = <P/^o 



(|)o sech^(KT)exp(ie,j,^), (12) 
i|/o sech^(KT)exp(i9y^), (13) 



where 9^ = 2k2 sgn(p20), 9^ + 2k2 (p2y/ |p20|) - 
ZQ (p2y — 2p2(j)) = are time independent constants. In our 
chosen frame of reference, these parametric soliton pulses 
have a group velocity of zero in the SVEA limit. 




FIG. 3; GFEA simulated simulton propagation: the pump field pro- 



file 

profile I 



for a 3-cycle pulse (A = 3) with a = 1. The signal field 
J)P is similar. 
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FIG. 4: GFEA simulated simulton velocity shifts — 5v as a function 
of bandwidth multiplier A. The pulses at A = 1 are approximately 10 
cycles long. The symbols represent the velocity shifts obtained by 
numerical simulation, the solid lines the corresponding theoretical 
prediction. The numerical labels give the group/phase velocity ratio 
a corresponding to the nearest data points. 



A. GFEA Simulations 



We use the same basic code as for the Kerr soliton simu- 
lations, with the different form of nonlinearity. We use the 
normalised units described after eqn.( 11). 

Fig. ^ shows a typical result for the propagation of a para- 
metric soliton. As for the Kerr soliton, few-cycle effects mod- 
ify the group velocity. Results for different soliton widths A 
and different a values are summarised in fig. Q The velocity 
shift remains remarkably linear even for extremely wideband 
pulses. Note that A = 5 corresponds to a 2 cycle pulse, for 
which the electric field profile would not appear particularly 
sech^ shaped, because of the small number of caiTier oscilla- 
tions. In reality, however, the parabolic form of the material 
dispersion assumed in the simulations will not be maintained 
over such a wide bandwidth, and other distortions are likely 
to predominate over few-cycle effects in these circumstances. 
In addition, as the bandwidth increases, the spectra of the two 
pulses will eventually overlap, despite the separation of their 
carrier frequencies. 



B. Theoretical Group Velocity Shift 

As described in section ^ we followed the method of 
Biswas and Aceves[|lO| to show how few-cycle corrections for 
NLS solitons would modify the group velocity. We now show 



that the same method can predict parametric soliton group ve- 
locity shifts. Eqn.(13) from Biswas and Aceves[p^ applied to 
the equations for parametric solitons gives the velocity shifts 



d J:^^^*^dT: 
dt, /+J(^*(|)t/T ' 
d J^^z\\t*\\t dz 



(14) 
(15) 



which we can evaluate by inserting the SVEA soliton pro- 
files (from eqn.(p^) and (|3])) along with the few-cycle pertur- 
bations to the propagation: 



Ri, = (2-o)ax\|/(j)7c0y, 



(2-o)8x^V(2cOv)' 



(16) 
(17) 



Solving to first order in the perturbation(s) and using the 
intermediate quantity ^ = J^^ld^^ i*^*^) di gives 



2 — o , —1671 
^ = ^:^(¥o(l)o(^o + v|/o(l)o(l)o)- 



V'lh = — - 



0)0 

8 2-a 
15 (Oif 



45 ' 



(¥o + ¥S): 



(18) 
(19) 
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and similarly using = J^Z'ld^ (¥*¥) dz, 

1 2-a /,,7 ,\ -1671 



IV. CONCLUSION 



2 COy 

4 2-0 



15 CO, 



45 



(20) 
(21) 



In our numerical simulations we used \\to — 3, <^o — —6i, 
and COy = 2cO||, = 20jt/A, so that 



5v = vy = v,^w -0.1019(2-0) A. 



(22) 



Because the velocity shift for each pulse of the pair mak- 
ing up the parametric soliton is the same (at least for first or- 
der GFEA corrections), the two pulses remain co-propagating, 
and the soliton survives. 

The predictions of eqn.(p^ are plotted as solid lines in fig. 
^ and the agreement with the numerical simulations is seen to 
be remarkably good, even for wideband pulses (e.g. A > 4). 
However, the effect of the higher order GFEA contributions, 
present in the simulations but not in eqn.(^2|), start to become 
visible above a « 1 .60, since near a — 2 the first order few- 
cycle correction does not dominate. 



We have investigated two types of soliton propagation be- 
yond the standard SVEA regime both theoretically and nu- 
merically. The most important result is that, according to 
the GFEA theory, soliton propagation remains robust in the 
few-cycle regime. This is obviously encouraging for proposed 
applications involving ultrashort solitons - although for such 
wideband pulses, there are other complications beyond just 
the few-cycle ones examined in this paper 



The major effect of the shortening pulses is a group veloc- 
ity shift, despite the fact that the perturbation term does not 
look like a a simple group-velocity term. It is likely that the 
few-cycle "phase twist" added to the propagation also affects 
the other properties of soliton pulses, e.g. collisions, which 
has obvious potential implications for soliton-based ultrafast 
optical logic gates. 
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